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We study fluctuations of conserved charges including baryon number, electric charge, and 

strangeness as well as the correlations among these conserved charges in the 2+1 flavor Polyakov- 

■ Nambu-Jona-Lasinio model at finite temperature. The calculated results are compared with those 

00 ' 

£ — . ■ obtained from recent lattice calculations performed with an improved staggered fermion action at 

in 
d 



two values of the lattice cutoff with almost physical up and down quark masses and a physical 



' value for the strange quark mass. We find that our calculated results are well consistent with those 

o ' 

obtained in lattice calculations except for some quantitative differences for fluctuations related 
with strange quarks. Our calculations indicate that there is a pronounced cusp in the ratio of the 
quartic to quadratic fluctuations of baryon number, i.e. Xa IXi '> a ^ the critical temperature during 
the phase transition, which confirms that Xa IxE ls a useiu l probe of the deconfinement and chiral 
phase transition. 
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I. INTRODUCTION 



QCD thermodynamics, for example the equation of state of the quark gluon plasma 
(QGP), phase transition of the chiral symmetry restoration, the deconfmement phase transi- 
tion and so on, has been a subject of intensive investigation in recent years. On the one hand, 
the deconfmed QGP are expected to be formed in ultrarelativistic heavy- ion collisions jl-8| 
(for example the current experiments at the Relativistic Heavy Ion Collider (RHIC) and the 
upcoming experiments at the Large Hadron Collider (LHC)) and in the interior of neutron 



stars 



12| ; On the other hand, studying the thermodynamical and hydrodynamical behav- 



iors of the QGP, especially the deconfmement and chiral phase transitions, is an elementary 
problem in strong interaction physics. 

Lattice QCD simulation is a principal approach to explore the properties of strongly 
interacting matter and its deconfmement and chiral phase transitions. In the past years, 
this method has provided us with lots of information about the QCD thermodynamics and 
phase transition from the confined hadronic phase to the deconfmed QGP one at finite 
temperature and limited chemical potential (see, for example, Ref. 13|-|29|). In response 
to the lattice QCD simulations, many effective models have been developed to describe 
the behavior and properties of strongly interacting matter, give physical interpretation of 
the available lattice data, and further to make predictions in the regions of phase diagram 
that can not be reached by the lattice calculations. The validity of these effective models 
is expected, since what governs the critical behavior of the QCD phase transition is the 
universality class of the chiral symmetry which is kept in these effective models. 

It has been known that fluctuations of conserved charges, for example the baryon num- 
ber, electric charge, and strangeness, are particularly sensitive to the structure and behavior 
of the thermal strongly interacting matter 30H32l| . Enhanced fluctuations are close related 



with the critical behavior of the QCD thermodynamics and phase transitions 33(. Further- 
more, the fluctuations and correlations of conserved charges and their high order cumulants 
provide information about the degrees of freedom (confined hadrons or deconfmed QGP) 



of strongly interacting matter at high temperature 



of the deconfmement and chiral phase evolution 2J-. 



34, 



351 ] . thus they are useful probes 
|. More important, the fluctu- 



ations and correlations of conserved charges can be extracted not only theoretically from 
the lattice QCD simulations and effective model calculations but also experimentally from 
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event-by-event fluctuations [32|, |35|, [37]. 



In this work, we will study the fluctuations and correlations of conserved charges and their 
high order cumulants in the 2+1 flavor Polyakov-Nambu-Jona-Lasinio (PNJL) model 38 ]. 
Compared with the conventional Nambu-Jona-Lasinio model, the PNJL model not only has 
the chiral symmetry as well as the dynamical breaking mechanism of thi s sy mmetry, but also 



include the effect of color confinement through the Polyakov loop [391447] . Our calculated 
results of the fluctuations and correlations of conserved charges will be compared with those 
obtained from very recent lattice calculations performed with an improved staggered fermion 
action at two values of the lattice cutoff with almost physical up and down quark masses 
and a physical value for the strange quark mass |28|]. On the one hand, our calculations 
provide insight about what can not be concluded from lattice calculations. For example, 
whether there is a cusp in the ratio of the quartic to quadratic fluctuations of baryon number 
Xa IX2i which is a valuable probe of the deconfinement and chiral dynamics, at the critical 
temperature can not yet be answered by lattice simulations, because of the larger errors in 
the calculations, but our calculations indeed find a pronounced cusp in this ratio xf/xf a t 
the critical temperature. Furthermore, our calculations give physical interpretation of the 
lattice data. On the other hand, comparing our calculated results with those obtained in 
lattice calculations, it allows us to test the validity of the 2+1 flavor PNJL model, and then 
provide the right direction for improving this effective model. 

The paper is organized as follows. In Sec. II we simply review the formalism of the 2+1 
flavor PNJL model. In Sec. Ill we introduce the fluctuations and correlations of conserved 
charges. In Sec. IV we give our calculated results of fluctuations of light quarks (up and 
down quarks) and strange quarks. In Sec. V we give our calculated results of fluctuations 
and correlations of conserved charges, and compare them with those obtained in lattice 
calculations. In Sec. VI we present our summary and conclusions. 



II. 2+1 FLAVOR PNJL MODEL 



In this work, we employ the 2+1 flavor Polya 
been discussed in details in our previous work ; 



cov-loop improved NJL model which has 
and the Lagrangian density for the 2+1 
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flavor PNJL model is given as 

8 

Cpnjl = 4> {ijfiD^ + 70/t - m ) ip + [{i> T a^) 2 + (V^TsTaVO' 

-K [det, (1 + 75) 1>) + det f (1 - T5 ) ^)) - U ($, $* , T) , (1) 

where ■?/> = {ip u ,'4 , d-i'4 , s) T is the three-flavor quark field, 

D" = & i - iA» with A" = ££A° , A = ^^°y = -iM- (2) 

A a are the Gell-Mann matrices in color space and the gauge coupling g is combined with 
the SU(3) gauge field A£(x) to define A^ix) for convenience. m = diag(mg, mf{, mg) is the 
three-flavor current quark mass matrix. Throughout this work, we take rriQ = itlq = m l , 
while keep being larger than m Q , which breaks the SU(3) / symmetry. The quark chemical 
potentials are contained in the matrix p, = diag(/x u , fid, /-0, and they can also be expressed in 
terms of chemical potentials for baryon number (/is), electric charge (/xq), and strangeness 
(Ms): 

12 11 ,11 

/x u = + -/i Q , = - and /x B = -/x B - -/x Q - /i S . (3) 

In the above PNJL Lagrangian, W ($,$*, T) is the Polyakov-loop effective potential, 
which is expressed in terms of the traced Polyakov-loop $ = (Tr c L)/N c and its conjugate 
$* = (Tr c L^)/iV c with the Polyakov-loop L being a matrix in color space given explicitly by 



L (x) = V exp 



f 

/ dr A A (x, t) 
Jo 



exp [if3A A ] , (4) 



with (3 = 1/T being the inverse of temperature and A<± = iA°. 

In our work, we use the Polyakov-loop effective potential which is a polynomial in $ and 



$* [421], given by 



$*$--(^ + $* d ) + -($*$) 2 , (5) 



T 4 2 6 v ' 4 

with 

T \ /To\ 2 , /T X 3 



6 2 (T) = a + ax 1^1 + a 2 ( + a 3 If . (6) 

Parameters in the effective potential are fitted to reproduce the thermodynamical behavior 
of the pure gauge QCD obtained from the lattice simulations, and their values are given in 
Table [H The parameter T is the critical temperature for the deconfinement phase transition 
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TABLE I: Parameters for the Polyakov-loop effective potential U 



a ai 02 03 b 3 64 



6.75 -1.95 2.625 -7.44 0.75 7.5 



to take place in the pure-gauge QCD and T is chosen to be 270 MeV according to the lattice 
calculations. 

In the mean field approximation, the thermodynamical potential density for the 2+1 
flavor quark system is given by 



p 



+ L l n [l _|_ 3$* e -(^|-M/)/T + 3$e -2(E/- M/ )/T + e -3(£/- W )/T 

3 L 

+ I. ^ [l + 3$ e -(£/+M/)/r + 3$ * e -2(B/+ W )/T + e -3(J5£+ W )/T 

3 

+2G(0 U 2 + 0/ + <p s 2 ) - AK<p u <p d S + $*, T). 



(7) 



where 



u, d, s) is the chiral condensate of quarks with flavor i, and we have the 



energy-momentum dispersion relation E % = \Jp 2 + Mf for its corresponding quasiparticle, 
with the constituent mass for the quark of flavor % being 



Mi = ml - AG(pi + 2K(j) 



•j <Pk- 



(8) 



Minimizing the thermodynamical potential in Eq. (J7j) with respective to <p u , (pd, <Ps, 3>, 
and $*, we obtain a set of equations of motion 



dtt on an 



9fi 



0. 



(9) 



c<^ u c<^ d dcj) s <9$ " <9$* 

Then, this set of equations can be solved as functions of temperature T and three flavor 
quark chemical potentials fi u , /i<j, and // s or chemical potentials of conserved charges fis, 
fi Q , and fis- 



III. FLUCTUATIONS AND CORRELATIONS 



Following the procedure in lattice calculations 



we focus on the derivatives of the 



pressure (P = — fl) of the thermodynamical system with respective to the chemical potentials 



corresponding to the conserved charge: baryon number, electric charge, and strangeness, i.e. 



BQS _ «9*+^(P/T 4 ) 



(10) 

Mfl,Q,S=0 



which are evaluated at Hb,q,s — 0- The x' s m Eq. ( TTUj) are in fact the generalized suscepti- 
bilities and they are related with the moments of charge fluctuations, i.e. SNx = Nx — (Nx) 
(X = B, Q, S) and the correlations among conserved charges. Taking the quadratic, quartic 
and the 6th order charge fluctuations for example, we have 

X? = ^Nx 2 ), (11) 
X? = ^r 3 ((SNx 4 )-3(6N x 2 ) 2 ), (12) 

,X _ 1 f/nr 6\ i c/rAr _A\/XAT. _2\ m«/\r 3\ 2 , 5n /i»r 2\3 



= ^(^(5iVx b )-15(5iV x 4 )(5iV^)-10(5iV^)^ + 30(5iV^) J j. (13) 
And the correlations among two conserved charges are 

Xfi y = ^(NxN Y ). (14) 

It should be noted that the fluctuations and correlations of conserved charges, i.e. the 
generalized susceptibilities in Eq. ( ITOj) are nonvanishing only when % + j + k is even at 
vanishing chemical potentials ^b,q,s — 0. 

In this work we use the method of Taylor expansion to compute the fluctuations and 
correlations of conserved charges in the PNJL model. Before our numerical calculations, we 
should determine the five parameters in the quark sector of the model. Values of the five 
parameters used usually in the literatures are obtained in Ref. [48j], which are m l = 5.5 MeV, 
m s = 140.7 MeV, GA 2 = 1.835, KA 5 = 12.36 and A = 602.3 MeV. They are fixed by fitting 
m n = 135.0 MeV, m K = 497.7 MeV, m n > = 957.8 MeV and f v = 92.4 MeV. For convenience 
to compare our calculations with the recent lattice simulations in which the strange quark 
mass has been tuned close to its physical value and the light quark masses have been chosen 
to be one tenth of the strange quark mass |28j, we employ the values of the five parameters 
given above except with m l = 14.0 MeV, which is consistent with the lattice calculations. 



IV. FLUCTUATIONS OF LIGHT QUARKS AND STRANGE QUARKS 



In order to understand the behavior of the fluctuations of conserved charges better, we 
will study the fluctuations of light quarks and strange quarks in this section. In Fig. [1] 
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FIG. 1: (color online). Quadratic (top), quartic (middle) and the 6th order (bottom) fluctuations 
of light quarks (here we denote them with \ u ) an d strange quarks as functions of the temperature 
(in unit of the critical temperature Tq) calculated in the PNJL model with m l = 14.0 MeV and 
ttiq = 140.7 MeV. We also show the results of the quadratic and quartic fluctuations obtained from 
calculations on 16 3 x 4 and 24 3 x 6 lattices in Ref. 28j|. Here SB denotes the corresponding limit 
value of Stefan-Boltzmann ideal gas. 

we show the quadratic, quartic and the 6th order fluctuations of light and strange quarks 
versus temperature. We find that the pseudo-critical temperature for the chiral restoration 
phase transition is 225 MeV and that for the deconfinement phase transition is 221 MeV 
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at vanishing chemical potential in the PNJL model with m l = 14.0 MeV, m s = 140.7 MeV 
and other parameters given above. More detailed discussions about the different phase tran- 
sitions can be found in Ref 



|. Since these two critical temperatures for the two different 
kinds of phase transitions are almost same, we just need one value of temperature to lo- 
cate the two phase transitions and here we choose Tq = 225 MeV which is used as unit in 
Fig. [TJ Furthermore, we also present the results of the quadratic and quartic fluctuations 



obtained from calculations on 16 3 x 4 and 24 3 x 6 lattices in Ref 28] in Fig. [TJ From the 
top panel of Fig. [1] ones can find that the quadratic fluctuations calculated in the PNJL 
model increase monotonously with the increase of the temperature, and the fluctuations 
of the heavier strange quarks are suppressed relative to those of the light quarks. These 



features are consistent with the lattice calculations (also the lattice simulations in Ref. |20|). 
Comparing our results with those from lattice calculations, we find that the quadratic fluc- 
tuations calculated in the PNJL model are relatively smaller than those calculated in the 
lattice simulations. And we also find that the quadratic fluctuations calculated in the lattice 
simulations rapidly approach the Stefan-Boltzmann ideal massless gas value, i.e. xf B = 1 
when the temperature is above the critical temperature, while the quadratic fluctuations 
calculated in the PNJL model approach the Stefan-Boltzmann value more slowly, which is 
because the current quark mass effect is obvious at high temperature in the PNJL model. 




FIG. 2: (color online). Ratio of s and u quark quadratic fluctuations as function of the temperature 



calculated in the PNJL model with m 



16 3 x 4 and 24 3 x 6 lattices in Ref. 



28]. 



14.0 MeV and mg 



140.7 MeV and in simulations on 



S 



In Fig. [2]we show the ratio of strange and light quark quadratic fluctuations as function of 
the temperature. Our calculations clearly indicate that the quadratic fluctuations of strange 
quarks are smaller than those of light quarks at low temperature, which are consistent with 
the lattice calculations. Furthermore, we find a plateau at the critical temperature in the 
curve of X2/X2 versus temperature in our calculations, as the solid line shows in Fig. |2J This 
plateau can also be found in the lattice data. 
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FIG. 3: (color online). Left panel: constituent masses of light quarks and strange quarks as 
functions of the temperature at vanishing chemical potentials in the PNJL model with m l = 
14.0 MeV and 7Bq = 140.7 MeV. Right panel: corresponding derivatives of constituent masses of 
light quarks and strange quarks with respective to the temperature, and here we use the absolute 
values of the derivatives. 



As for the quartic fluctuations, we find that the fluctuations of light quarks calculated in 
the PNJL model is consistent with those obtained from the lattice calculations, especially 
during the phase transition where a sharp peak appears at Tc- Furthermore, lattice calcula- 
tions show that there is a relatively lower peak at Tc as well for the quartic fluctuations of 
strange quarks, while our calculations indicate that the peak corresponding to the strange 
quarks is shifted to about 1.2Tq and becomes broader and lower. To explore the reason 
for the difference between our and lattice calculations, we show the constituent masses of 
light quarks and strange quarks, and also their derivatives with respective to temperature 
as functions of the temperature in Fig. [3j One can find that the constituent mass of light 
quarks rapidly decreases to their current quark mass at Tc, which corresponds to a sharp 
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peak of the curve \dM u /dT\ at the critical temperature as the right panel of Fig. [3] shows. 
However, for the strange quarks, the constituent mass still decreases slowly beyond the crit- 
ical temperature, which results in a plateau of the curve \dM s /dT\ at about 1.2Tc- That 
is the reason why the quartic fluctuations of strange quarks present a broad peak at 1.2Tq 
in the PN JL model. In the bottom panel of Fig. [TJ we show the 6th order fluctuations of 
light quarks and strange quarks. One can find that the 6th order fluctuations of light quarks 
oscillate violently at Tc, while those of strange quarks oscillate mildly at relatively larger 
temperature. 
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FIG. 4: (color online). Ratio of the quartic to quadratic fluctuations for light quarks (left panel) 
and strange quarks (right panel) as functions of the temperature calculated in the PNJL model 



with m\ 



14.0 MeV and mg 



140.7 MeV, which are also compared with those obtained from 



calculations on 16 3 x 4 and 24 3 x 6 lattices in Ref. 



28|]. 



In the following we pay more attentions to the ratio of the quartic to quadratic fluctuations 
of the quark number, since this ratio is believed to be a valuable probe of the deconfinement 
and chiral phase transitions 24j-|26l. |36| . Fig. H] shows this ratio for the light quarks and 
strange quarks calculated in our PNJL model, and we also show the corresponding results 
obtained in the lattice simulations for comparisons. Our calculations indicate that, for both 
light quarks and strange quarks, the ratio of the quartic to quadratic fluctuations of the 
quark number increases with the temperature below the critical temperature, while decreases 
when the temperature is above the critical one. So there is a peak in the curve of the ratio 
versus temperature at the critical temperature for both light quarks and strange quarks. 
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These features are consistent with those obtained in the lattice calculations. Comparing the 
two panels of Fig. H] we find that the consistency between our calculations and the lattice 
simulations for the light quarks is better than that for the strange quarks, which can be 
expected, since for the quartic fluctuations of strange quarks the difference between the 
PNJL model and the lattice simulations is quite larger than that for the light quarks, as 
shown in the middle panel of Fig. [TJ 

Next, we discuss the ratio of the quartic to quadratic fluctuations when the temperature 
is away from the phase transition temperature. As the temperature is high, especially when 
the temperature is quite above the critical temperature, the system can be approximated 
as noninteracting massless gases. That is the Stefan-Boltzmann limit, and in this limit the 
pressure of the quarks and antiquarks can be easily obtained as 



P 



f=u,d,s 



7tt 2 l//y\ 2 1 /HfY 
180 6 V T / 12n 2 \T J 



(15) 



So in the Stefan-Boltzmann limit, Xi/X2 = 6/7T 2 . 

1.2, . , . , , , . 1 10 




FIG. 5: (color online). Left panel: Polyakov-loop versus temperature in two cases at vanishing 
chemical potentials in the PNJL model with m l = 14.0 MeV and mg = 140.7 MeV ($* = <& 
at vanishing chemical potentials). The blue dashed line corresponds to our full calculations, and 
the red solid line to the calculations in which and <&* are set to be zero artificially when the 
temperature is below the critical temperature. Right panel: ratio of the quartic to quadratic 
fluctuations for light quarks corresponding to the two cases. Therefore, the blue dashed line 
represents the same results which have been shown in the left panel of Fig. 
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On the contrary, when the temperature is below the critical temperature, the constituent 
masses of quarks are much larger than their current masses as the left panel of Fig. [3] shows. 
Furthermore, quarks tend to be combined to form colorless hadrons, which is the color 
confinement. In the PNJL model, we can not include the degrees of freedom of hadrons, but 
this effective model embodies the confinement effect by suppressing the excitations of the 



one and two quark states and only permitting those of three quark states 38]. In details, 
when the temperature is much lower than the critical temperature, the Polyakov-loop and its 
conjugate $* — > 0, and then most contributions to the pressure of the thermodynamical 
system come from the three quark states. Since the constituent masses of quarks are much 



larger than the temperature, we can use the Boltzmann approximation 36), and the pressure 
is 



f=u,d,s 

Where Ki is a Bessel function. Then we can easily obtain Xi/X2 = 9. It should be em- 
phasized that, since the temperature driven deconfinement transition is not a strict phase 
transition but a continuous crossover 2l|, |22j], the Polyakov-loop $ (or $*) is not vanishing 



when the temperature is below the critical one as the blue dashed line shows in the left 
panel of Fig. |5j Assuming $ = $* = 1 even when T < Tq, we obtain 

f=u,d,s 

Then we have x^lxi = 1 in this situation. Since $ (or $*) is between and 1 at low 
temperature as we have mentioned above, we expect x^lxi to be between 1 and 9 when 
T < Tc- This expectation is verified in our numerical calculations as shown in Fig. HI 
For comparison, we set $ and $* to be zero artificially when the temperature is below the 
critical temperature as the red solid line shows in the left panel of Fig. then we calculate 
the corresponding Xa/X2- The results are shown in the right panel of Fig. [51 which clearly 
indicates that X4/X2 approaches 9 when T < Tq. One can easily find that the ratio of the 
quartic to quadratic fluctuations calculated in the case in which the Polyakov-loop is set to 
be zero artificially is not consistent with the lattice calculations. 

So far, we have found that the calculations in the PNJL model for the light quark 
fluctuations are consistent with those in the lattice simulations, while for the strange 
quark fluctuations, since the current mass of the strange quarks is large, there are rela- 
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FIG. 6: (color online). Left panel: quartic fluctuations of light quarks and strange quarks as 
functions of the temperature in the PNJL model with different quark current masses, one with 



14.0 MeV and ttt-q = 140.7 MeV (same as our previous calculations) and the other with 
3.79 MeV and rrig = 104 MeV. In the same way , the corresponding results obtained in QCD 

281 ] are also presented. Right panel: ratio of the 



simulations on 16 3 x 4 and 24 3 x 6 lattices in Ref. 
quartic to quadratic fluctuations for strange quarks as functions of the temperature in the PNJL 
model with the two sets of current quark masses. 

tively larger differences between the results in the PNJL model and those in lattice cal- 
culations. Considering the value of the strange quark current mass used in the PNJL 
model (wq = 140.7 MeV) is relatively larger than the middle value 104 MeV given by 



the PDG2008 {ml 



1041^6 MeV 



49|), we set m s = 104 MeV and m l Q = 3.79 MeV 
in the PNJL model and leave other parameters unchanged (here m l = 3.79 MeV is 
the middle value of the average mass of up and down quarks given by the PDG2008 
(m = (m u + md)/2 = 3.79^Jg MeV 49J)). Then we repeat the calculations above. Fig. [6] 
shows the quartic fluctuations of light quarks and strange quarks, and the ratio of the quar- 
tic to quadratic fluctuations for strange quarks with the two sets of current quark masses. 
We find that the height of the peak in the quartic fluctuations of light quarks at the critical 
temperature is enhanced due to the reduction of light quark current mass, and with the 
decrease of the strange quark mass, the broad peak of the curve xt move a little to low tem- 
perature. However, there are also differences between the PNJL model with m s Q = 104 MeV 
and the lattice simulations for the xt an d X s Jx\- 
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V. FLUCTUATIONS AND CORRELATIONS OF CONSERVED CHARGES 




FIG. 7: (color online). Quadratic (top), quartic (middle) fluctuations of baryon number, and their 
ratio (bottom) as functions of the temperature calculated in the PNJL model with m l = 14.0 MeV 
and rag = 140.7 MeV. The corresponding results obtained from calculations on 16 3 x 4 and 24 3 x 6 



lattices in Ref. 



281] are also presented. 



We have discussed the fluctuations of light quarks and strange quarks in details above. 
In the following, we will focus on the fluctuations of baryon number, electric charge, and 
strangeness. 

Fig. [7] shows the quadratic and quartic fluctuations of the baryon number, and their ratio 
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versus temperature calculated in the PNJL model and in the lattice simulations. Same as the 
fluctuations of light quarks and strange quarks, We find that the quadratic fluctuations of 
baryon number increase monotonically with the temperature, while the quartic fluctuations 
develop a cusp at the critical temperature and decrease with increase of the temperature 
when T > Tq- As a matter of fact, the singular behavior of the baryon number fluctuations 
is expected to be controlled by the universal 0(4) symmetry group at vanishing chemical 



potential and vanishing light quark mass 31] . And the baryon number fluctuations are 
expected to scale like 



T — Tn 2-n-a 

(18) 

J-c 

with a ~ —0.25. As Eq. (1181) shows, with the increase of the order, the singular properties 
of the baryon number fluctuations become more and more prominent, which are confirmed 
both in the calculations of the PNJL model and in the lattice simulations. Comparing the 
results obtained in the PNJL model and those in lattice calculations, we find that our results 
are consistent with the lattice results qualitatively, but the differences in quantities between 
the two results still exist, especially at high temperature. The reason for these differences 
has been mentioned above, i.e. when the temperature is high, the current quark mass effect 
becomes more and more important. 

The bottom panel of Fig. [7] shows the ratio of the quartic to quadratic fluctuations of 
the baryon number. We find that this ratio approaches 1 at low temperature, which is 
consistent with the lattice results as well as the calculations in the hadron resonance gas 
model 
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50| . An interesting thing is that we find a pronounced cusp in the ratio xf / X2 a ^ 
the critical temperature in our calculations. However, whether the cusp exists in the lattice 
results is not clear, because the errors are very large at low temperature as the figure shows. 
We have emphasized that answering the question whether there is a cusp in the ratio of the 
quartic to quadratic fluctuations is very important, because this ratio is a valuable probe of 
the deconfinement and chiral phase transitions and have received lots of attentions recent 



years 
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361 ] . Our calculations indicate that the cusp in the ratio X?/xf is prominent even 
with physical strange quark mass rriQ = 140.7 MeV and relatively larger (comparing with 
the physical value) light quark mass m l = 14.0 MeV. 

In Fig. El we show the quadratic, quartic and the 6th order fluctuations of electric charge 
versus temperature obtained in our calculations and lattice simulations. Once more, we find 
that the electric charge quadratic fluctuations increase monotonically with the temperature 

15 




FIG. 8: (color online). Quadratic (top), quartic (middle) and the 6th order (bottom) fluctuations 
of electric charge as functions of the temperature obtained in the PNJL model with rn l = 14.0 MeV 
and uiq = 140.7 MeV and in the QCD simulations on 16 3 x 4 and 24 3 x 6 lattices in Ref. [28.]. 



and the quartic fluctuations of the electric charge present a prominent cusp at the critical 
temperature. As for the 6th order fluctuations of electric charge, we find a pronounced 
oscillation happening during the deconfmement and chiral phase transitions. x& * s vanishing 
when the temperature is below about 0.8T C , and with the temperature being increased 
beyond 0.8Tc the 6th order fluctuations increase rapidly. When the temperature is near 
Tc, Xe changes from a positive value to a negative one abruptly. Then, Xq increases with 
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the temperature and approaches the Stefan-Boltzmann value (xq' SB — 0) from below. From 
the bottom panel of Fig. [BJ One can find that the 6th order fluctuations of electric charge 
are almost odd functions with respective to T — Tc- Comparing Fig. [HJwith Fig. we find 
that the fluctuations of electric charge calculated in the PNJL model are better consistent 
with those obtained in lattice calculations than the fluctuations of the baryon number. This 
is because more contributions to the electric charge fluctuations come from the light quarks, 
since the absolute value of the electric charge of u quarks 2/3 is larger than that of s quarks 
1/3. 
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FIG. 9: (color online). Left panel: Ratio of the quartic to quadratic electric charge fluctuations 
versus temperature obtained in the PNJL model with wL = 14.0 MeV and = 140.7 MeV and 



in the QCD simulations on 16 3 x 4 and 24 3 x 6 lattices in Ref 



281 ] . Right panel: Ratio of the 6th 



to 2nd order fluctuations of electric charge versus temperature calculated in the PNJL model with 



14.0 MeV and = 140.7 MeV and in the QCD simulations on 16 3 x 4 lattices in Ref. 



28]. 



In Fig. [9j we show the ratios of the 4th to 2nd order and the 6th to 2nd order fluctuations 
of electric charge as functions of the temperature obtained in the PNJL model and the lattice 
simulations. One can find that the results given by the PNJL model are well consistent with 
those obtained in the lattice simulations. Comparing X4/X2 with X4/X2 m the bottom 
panel of Fig. [7J we find there is a prominent cusp at the critical temperature in the curve 
of xi IX2 as same as that in the curve of xf /xf • However, different from xf/xf > a t l° w 
temperature X4/X2 * s larger than unity and is a raising function of the temperature, since 
it is expected that at low temperature quartic fluctuations of electric charge are enhanced 
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relative to the quadratic fluctuations [28J. As for the ratio Xs/x-zi we find it oscillates 
violently at the phase transition temperature. 
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FIG. 10: (color online). The 6th order fluctuations of strangeness (left panel) and the ratio xf/xf 
(right panel) as functions of the temperature calculated in the PNJL model with m l = 14.0 MeV 
and rag = 140.7 MeV and in the QCD simulations on 16 3 x 4 lattices in Ref. 28^ . 



We have discussed the quadratic and quartic fluctuations of the strange quarks in details 
above. It is easily verified that the fluctuations of strange quarks are identical to those of 
strangeness. Therefore, we will focus on the 6th order fluctuations of strangeness in the 
following. Fig. [10] shows the 6th order fluctuations of strangeness and the ratio xf / xf as 
functions of the temperature calculated in the PNJL model, which are also compared with 
the QCD simulations on 16 3 x 4 lattices in Ref [28(. We find that xf an d xf/xf obtained 
in the PNJL model are qualitatively consistent with those obtained in lattice simulations , 
i.e. xf an d xf/xf oscillate during the deconfinement and chiral phase transitions. However, 
there are some quantitative differences between these two approaches: the temperature 
where xf ano ^ X&lxt vanish during the oscillation in the PNJL model is shifted from Tq 
to about 1.2Tc with respective to the results in the lattice simulations. Similar phenomena 
have been found in the calculations of xl an d X4/X2 above in Fig. [61 since much larger 
current mass of the strange quarks results in that the chiral restoration phase transition for 
the strange quark sector occurs at relative larger temperature than that for the light quark 
sector as Fig. [3] shows. 

In Fig. [TTJ we show correlations among conserved charges, in more details, including 
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FIG. 11: (color online). Correlations of baryon number and electric charge normalized to quadratic 
fluctuations of baryon number (top), baryon number and strangeness to quadratic fluctuations of 
baryon number (middle), and electric charge and strangeness to quadratic fluctuations of electric 
charge (bottom) as functions of the temperature obtained in the PNJL model with m l = 14.0 MeV 
and ?7t,q = 140.7 MeV and in the QCD simulations on 16 3 x 4 and 24 3 x 6 lattices in Ref. [201 . 



the correlations of baryon number and electric charge, baryon number and strangeness and 
those of electric charge and strangeness. The former two correlations are normalized to the 
quadratic fluctuations of baryon number and the last one is to quadratic fluctuations of 
electric charge. One can find that Xn /xf calculated in the PNJL model decreases with the 
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increase of the temperature, which agrees with the calculations in lattice simulations. When 
the temperature is about 1.5Tq, this ratio approaches zero, the limit value of noninteracting 
massless quark gas, since the sum of electric charges of light quarks and strange quarks 
vanishes in this limit. Ones can also find that our calculated results of X11/X2 are a 
little larger than those of lattice calculations at the critical temperature, but well consistent 
with the lattice results at high temperature. The correlation between baryon number and 
strangeness has been thought as a useful diagnostic of strongly interacting matter, which 
can be extracted not only theoretically from lattice QCD and effective model calculations 
but also experimentally from event-by-event fluctuations [35(]. We find that the correlation 
of baryon number and strangeness approaches zero at low temperature, since the strange 
hadrons (here in the PNJL model are the three quark clusters with strange quarks) are 
much heavier to excite in the thermodynamical system. The ratio —X11/X2 calculated in 
the PNJL model increases with the temperature and approaches the Stefan-Boltzmann value 
at about l.bTc, which agrees with the lattice QCD simulations. Same as the correlation of 
baryon number and strangeness, the correlation of strangeness and electric charge is also 
suppressed at low temperature and the ratio Xii /xi increase with the temperature. One 
can find that our calculated results of X11 /X2 are wen consistent with those of the lattice 
simulations. An interesting thing is that both the results given by the two approaches 
indicate that there is a plateau in the ratio of X11 /X2 a ^ ^ e critical temperature. 



VI. SUMMARY AND DISCUSSIONS 



In this work, we have studied the fluctuations and correlations of the conserved charges, 
i.e. the baryon number, electric charge and the strangeness, in the Polyakov-Nambu-Jona- 
Lasinio (PNJL) model at finite temperature. We have also compared our calculated results 
with those obtained from the recent lattice calculations performed with an improved stag- 
gered fermion action at two values of the lattice cutoff with almost physical up and down 
quark masses and a physical value for the strange quark mass 28]. 

In order to study the fluctuations of conserved charges much better, we have calculated the 
fluctuations of light quarks (up and down quarks) and strange quarks before the calculations 
for the conserved charges. Our calculations indicate that the quadratic fluctuations of light 
quarks, strange quarks and conserved charges (including the baryon number, electric charge 
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and the strangeness) increase from zero at low temperature to their corresponding Stefan- 
Boltzmann values at high temperature, since particles excited at finite temperature are 
increased with the increase of temperature; the quartic fluctuations are characterized by an 
pronounced cusp (the cusp is high and sharp for light quarks and low and broad for strange 
quarks, but both are prominent) during the deconfinement and chiral phase transitions; the 
6th order fluctuations for the quarks and conserved charges oscillate with the occurrence of 
the phase transitions. These qualitative features shown in our calculations are well consistent 
with those given by the lattice calculations, which confirms that the effective model (here is 
the PNJL model) captures the right symmetry of the QCD, i.e. the chiral symmetry group, 
since different phase transitions are govern by different symmetry groups. 

Comparing our calculated results with those obtained in lattice QCD simulations quan- 
titatively, we have found that the fluctuations of light quarks obtained in the PNJL model 
are well consistent with those calculated in the lattice simulations. However, for the fluc- 
tuations of strange quarks, there are discrepancies between these two approaches. In more 
details, calculations of the PNJL model indicate that the cusp of the quartic fluctuations 
and the oscillation of the 6th order fluctuations for strange quarks are located at about 
1.2Tc, while they are still located at Tq as same as the light quarks in the lattice calcula- 
tions in Ref. 28j. The reason for this discrepancy is that the pseudo-critical temperature 



of the chiral restoration phase transition for strange quarks is a little larger than that for 
light quarks in the PNJL model due to the larger current quark mass of strange quarks 



and locations of the cusp of the quartic fluctuations and the oscillation of the 6th order 
fluctuations for strange quarks is related with this pseudo-critical temperature for strange 
quarks, which is about \flTg in the PNJL model. On the contrary, results obtained in 
lattice calculations in Ref. [28| indicate that the pseudo-critical temperature of the strange 
quark chiral phase transition is same as that for light quarks. However, we should emphasize 
that in some other lattice simulations, for example in Ref. 22], these two pseudo-critical 
temperatures are different and the pseudo-critical temperature for strange quarks is larger 
than that for light quarks. Furthermore, because of the effects of current quark masses 
in the PNJL model, especially for the strange quarks, fluctuations obtained in the PNJL 
model approaches to their corresponding Stefan-Boltzmann values more slowly than those 
obtained in lattice calculations. Comparing the fluctuations of electric charge with those of 
baryon number, we find that the former calculated in the PNJL model are better consistent 
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with those obtained in lattice calculations than the latter, since the fluctuations of electric 
charge have more contributions from light quarks than the fluctuations of baryon number. 
Although the conclusion about whether there is a cusp during the phase transition in the 
ratio of quartic to quadratic fluctuations of the baryon number, i.e. X4/X2 > ^ s n °t c l ear i n 
lattice calculations due to the large errors in their simulations, We indeed find a pronounced 
cusp in this ratio Xa / xf a ^ the critical temperature in our calculations, which confirms that 
xf 7X2 * s a valuable probe of the deconfinement and chiral phase transitions. 

Furthermore, we have also studied the correlations among conserved charges in the PNJL 
model and compared our calculated results with those obtained in lattice simulations. Same 
as the fluctuations of conserved charges, except for some small quantitative differences, these 
correlations among conserved charges calculated in the PNJL model are well consistent with 
those obtained in the lattice calculations. 
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